Synchronization of Coupled Systems with Spatiotemporal Chaos 
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We argue that the synchronization transition of stochastically coupled cellular automata, dis- 
covered recently by L.G. Morelli et al. (Phys. Rev. 58 E, R8 (1998)), is generically in the directed 
percolation universality class. In particular this holds numerically for the specific example studied by 
these authors, in contrast to their claim. For real- valued systems with spatiotemporal chaos such as 
coupled map lattices, we claim that the synchronization transition is generically in the universality 
class of the Kardar-Parisi-Zhang equation with a nonlinear growth limiting term. 

PACS numbers: 05.70.Jk, 61.25.Hq 



Since the pioneering work of Fujisaka et al (Q-Q and 
others f|-Q], synchronization of chaotic systems has be- 
come a very intensely studied subject, partially due to 
hopes that this could lead to applications in control and 
secure communications ||. 

In spatially extended systems, synchronization can 
appear in (at least) two forms. On the one hand, one can 
ask whether distant regions in a single such systems can 
oscillate in phase. After this phenomenon was observed 
by Chate et al || in high dimensional cellular automata, 
it was realized that it can be mapped onto the Kardar- 
Parisi-Zhang (KPZ) problem of the growth of a random 
surface [ [To] , with a synchronized state corresponding to 
a globally smooth phase <j) = 0(x) [|ll] p"3[ . 

In the present paper we shall deal with another prob- 
lem, namely that of mutual synchronization of two iden- 
tical locally coupled systems. If these systems are de- 
scribed by state variables x\ and y\ (for simplicity we 
assume here discrete 1-dimensional space i and discrete 
time t), we write the evolution in general as 
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The function / is nonlinear such that the evolution is 
chaotic for e = 0. Due to sensitive dependence on initial 
conditions, x* and y* will be completely uncorrelated in 
this case, unless they started with identical initial condi- 
tions. Synchronization should only be expected for e > 
if g(x) is negative for positive small x, so that any small 
difference x\ — y\ will be damped by the last terms in 
eq.(l). 

For chaotic systems with a finite number of degrees 
of freedom, there is a finite synchronization threshold 
e c , with intermittent behavior and 'riddled' fl4|| attrac- 
tor basins near e — e c [0 HI. Recently, it was found 
numerically that essentially the same phenomena occur 
in spatially extended systems jl^-|l|. While chaos in 
systems with few degrees of freedom requires x to be 
real-valued, spatio-temporal chaos can occur also in sys- 
tems with discrete x\, co-called cellular automata (CA). 



Therefore one can ask whether the phenomenon of mu- 
tual synchronization can occur also in the latter p5[ , and 
whether there are universal scaling laws at the synchro- 
nization threshold which apply both to real-valued sys- 
tems (coupled map lattices and partial differential equa- 
tions) and to CA. While the former question was asserted 
positively in |f7j , we shall argue that the latter has a neg- 
ative answer. For CA, the synchronization threshold is 
generically in the directed percolation universality class 
p0| , while it is for continuous systems in the universal- 
ity class of KPZ growth with a nonlinear growth limiting 
term 

Let us first study the case of 1-d cellular automata. 
The specific system studied in [[L7[ was two copies evolv- 
ing according to Wolfram's p4| rule 18 with periodic 
boundary conditions, and endowed with an additional 
stochastic coupling term. In rule 18, x\ can assume 
two values or 1, and the evolution function / de- 
pends only on x\ itself and its two nearest neighbors, 
/(0,0,1) = /(1,0,0) = 1 and f(x i - 1 ,x i ,x i+1 ) = else. 
The coupling was realized as follows: after applying the 
above rule to both x and y, it was checked whether 
Xi = Hi- If not, a random number is drawn uniformly 
from [0, 1]. If this number is less than some fixed num- 
ber p, a second random number is drawn and, depend- 
ing on that, either Xi is put equal to j/j or g/j is put 
equal to Xi. Thus Xi — yi is enforced with probability 
p, while both Xi and yi are left untouched with proba- 
bility 1 — p. It was found numerically that the system 
synchronizes for p > p c — 0.193. For p < p c the density 
of sites with Xi =^ g/j scales for t — > oo as (jp c — p)@ with 
(3 = 0.34, while it decays for p > p c with a characteristic 
time T which scales as T ~ (p — p c ) _!/ " with v\\ = 1. 
Since these exponents disagree grossly with the DP val- 
ues (3 = 0.2765, U\i = 1.7338 Q, it was concluded that 
this transition is not in the DP universality class. 

Unfortunately, the above estimates are flawed for sev- 
eral reasons. The first is that it is notoriously difficult to 
measure f3 directly in DP and similar processes, due to 
the very slow convergence towards the stationary state. 
At p = p c , the density of 'active' sites in DP (correspond- 
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ing to sites with Xi ^ yi in the present model) scales as 
p ~ t~ s with S = 0.1595 [p5[ . Near p c , this means that 
one has to wait excessively long until the stationary state 
is reached. Much more reliable results are obtained by 
following the approach towards the asymptotic state, e.q. 
by measuring the decay of p with time. 

The second problem is that rule 18 is well known to 
have very slow convergence towards its asymptotic state 
pfjj, in contrast to claims made in pjj]. Therefore the 
strategy of discarding a transient of a few hundred time 
steps used in |17j] is bound to induce errors. When start- 
ing with a random initial state, rule 18 orders into do- 
mains in which x\ is zero either for even i + t or for odd 
i + t. The boundaries between these domains move ac- 
cording to annihilating random walks, so that the domain 
sizes grow ~ \ft and the density of domain walls decreases 
as Asymptotically, the entire lattice is one single 

domain. On the sublattice where x\ is not identically 
zero, its evolution follows the 'additive' rule 90 given by 
/(0,0,1) = /(1,0,0) = /(0,1,1) = /(1,1,0) = 1. The 
invariant state of the latter is completely random. There- 
fore rule 90 must show the same synchronization thresh- 
old as rule 18 and the same critical exponents, but it 
involves no transient whatsoever if one starts with ran- 
dom initial conditions. 
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FIG. 1. Log-log plot of p(t) for rule 90, for several values of 
p: 0.1902, 0.19059, 0.19065, 0.1910 from top to bottom. Sta- 
tistical errors are smaller than the thickness of the lines. The 
straight dashed line has slope —0.1595 as predicted by DP. 



The density p{t) for rule 90 is shown in Fig.l for sev- 
eral values of p. To obtain these data we used ca. 1000 
lattices of size L = 10000 for each p, which gav e a sam- 
ple more than 100 times larger that that of |ft7|] . We see 
clearly a power behavior for large t (with strong small- 
t corrections) for p = p c = 0.19061 ± 0.00003. This is 
quite far from the value given in fl7j| and implies im- 
mediately that the inserts in Figs. 2 and 3 of that pa- 
per are very misleading. We also see from Fig.l that 



S = 0.159 ± 0.003 in excellent agreement with DP. Af- 
ter having determined p c in this way, we performed very 
long runs (with up to 400,000 time steps and with L up 
to 40,000) for p < p c in order to estimate p(t = oo). Such 
long runs were needed since otherwise we would have suf- 
fered from systematic errors. Results are shown in Fig. 2 
and give (3 — 0.277 ± 0.007, again in perfect agreement 
with DP 

Finally, we show in Fig. 3 the density decay p(t) for 
rule 18, with random initial conditions and without dis- 
carding any transient. We see indeed rather large correc- 
tions to scaling for times ^ 10 2 . If we would locate the 
critical point by means of least square fits including the 
short time region, we would systematically overestimate 
p c and 5, just as was done in |]l7|. 



0.6 
0.5 



0.4 



0.3 



0.2 



0.1 




0.001 



0.01 



Pc-P 



FIG. 2. Log-log plot of p(t = oo) against p c — p, with p c 
as obtained from Fig.l. The dashed line has slope 0.2765 as 
predicted by DP. 
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FIG. 3. Same as Fig.l, but for rule 18 and for two values 
of p only: 0.1906 (top) and 0.192 (bottom). The first gives 
perfect agreement with DP for large t but large deviations 
at small t. The second would give a better least square fit to 
a straight line for 10 < t < 5000, but this would yield wrong 
estimates of p c and S. 
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We also studied rule 22. This is not an additive 
rule, and it has a nontrivial invariant measure with 
maybe zero entropy, but non-zero Lyapunov exponent 
PH (for the notions of chaos, entropy and Lyapunov ex- 
ponents for CA see [^9| ) . Even if the entropy is not zero, 
there are very long ranged correlations in the invariant 
measure of rule 22 [ pSfl . It is thus of interest to see 
whether there is still a synchronization transition, and 
whether it is still in the DP universality class. We found 
again perfect agreement with DP. The critical point is at 
p c = 0.22735 ± 0.00005. It is easily seen that the sum 
of right- and left-moving Lyaponov exponents has to be 
positive for p c to be non-zero in any 1-dimensional CA. 
But the above values for p c and the known values for the 
Lyaponov exponents for rules 22 and 90 J||,^j| suggest 
that there is no simple relationship beyond this qualita- 
tive criterion. 

On the theoretical side we also have no formal proof of 
the universality with DP, but we can use exactly the same 
heuristic arguments which were used in [ |30[ to argue that 
damage spreading transitions are generically in the DP 
universality class. We refer to [|o| for a detailed discus- 
sion, including caveats and limitations of the expected 
universality. Our present results underline again the re- 
markable robustness of DP critical behavior. In contrast 
to a statement made in |i~7| , up to now DP universality 
was verified in all tested cases (even if the original au- 
thors often found violations, such as in the present case), 
provided the criteria listed in p0[ were met. 

Let us now shortly discuss systems with continuous 
variables such as coupled map lattices. The main dif- 
ference between these and CA is that synchronization is 
never perfect for finite time, even if e > e c . Instead, the 
differences \x\ — y\\ decrease exponentially with t when 
the systems synchronize. But this means that close to 
threshold statistical or chaotic fluctuations can make the 
system de-synchronize again, at least locally. Technically 
spoken, the system does not enter an absorbing state 
when it synchronizes locally, in contrast to the discrete 
case. This implies very different scaling exponents, as 
first noticed in [£l| and verified in |o]|2|j2|]. The generic 
stochastic partial differential equation with these features 
contains a diffusion term, a local nonlinear term, and a 
multiplicative noise term. The logarithm of the field ap- 
pearing in this stochastic PDE satisfies the KPZ equa- 
tion with an additional term which prevents the height 
variable from overcoming a barrier which we can conve- 
niently place at h = 0. The synchronization transition 
in this version corresponds to a transition from a sur- 
face pinned at h ss (desynchronized state), to a surface 
drifting towards h = — oo (synchronized state). 

We conjecture that the transition found in the neu- 
ral network model of [|l8| is in this universality class. It 
would be interesting to make detailed simulations of that 
model to verify this numerically. 



On the other hand we conjecture that the synchro- 
nization transition studied in Jigj] is not an this univer- 
sality class. In [^9) the coupling strength was called 7, 
and perfect synchronization occurred exactly at 7 = 1. 
As seem from eq.(l) of ]l9| ], the equivariance group of the 
coupled system changes at 7 = 1. For 7^1 the system is 
invariant under phase transformations A1.2 ~ > A\^ 
and under the exchange A\ <-> A 2 . For 7 = 1 one has the 
additional symmetry under phase rotations A\ ± A 2 — > 
(Ai ± A 2 )e t ' t ' ± . As stressed in (de-)synchronization 
is essentially a phenomenon of spontaneous symmetry 
breakdown. As in other transitions with spontaneous 
symmetry breaking, the universality class of such tran- 
sitions should depend crucially on the type of symmetry 
broken, and should be particular sensitive to symmetry 
changes at the critical point. 
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